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Abstract 

The dominant contribution to the inclusive decays of B mesons into the charm- 
strangeness mesons Dg and D* is expected to be given by the partonic process 
6— >c+ {D^ , D*^). We determine the nonperturbative 0(l/m^) and the 0{as) 
radiative corrections to 6— > c + {D~,Dl~) and thereby to the inclusive decays 
B — > Xc + {Dj , Dl~). The new feature of our calculation is that we separately 
determine the nonperturbative and the 0{as) corrections to the longitudinal (L) 
and transverse (T) pieces of the spin 1 D*~ meson. The longitudinal/transverse 
composition of the D*~ can be probed through its two principal decay modes 
D*~ + 7 and Dl~ -^D~ + vr" for which we write down the angular decay 

distributions. 



1 Introduction 



In a recent paper Aleksan et al. have convincingly argued that the inclusive decay 
B Xc + {Dj,D*~) is dominated by the partonic process 6 — >• c + {Dj,D*~) |^. 
The basic assumption is that factorization holds for the nonleptonic decay process 
B ^ Xc + {D~,D*~). One can then factorize the transition into a current-induced 
B ^ Xc transition and a current-induced vacuum one-meson transition. The leading 
order contribution to the B ^ X^ transition is given by the partonic 6 — >■ c transition. 
Corrections to the leading order result set in only at 0(l/m^). They can be estimated 
using the methods of the operator product expansion in HQET. 

Aleksan et al. also pointed out that it would be interesting to experimentally mea- 
sure the longitudinal/transverse composition of the spin 1 meson D*~ in this inclusive 
decay which they computed at the Born term level. In an accompanying paper the 
same authors calculated the 0{as) corrections to the inclusive rates into the spin 
D~ and the spin D*^ without ,[| however, separating the longitudinal (L) and 
transverse (T) contributions in the spin 1 case. It is the purpose of this paper to fill 
the gap left by and to provide analytical formulae as well as the relevant numerical 
results for the 0{as) L/T content of the D*~ in this reaction. We emphasize that 
the 0{as) corrections calculated here and in are only partial. There are also non- 
factorizing 0{as) corrections as Beneke et al. [^] and Chay have explicitly shown for 
the exclusive decays B—^ini and B^D'^*^n^, respectively. The non-factorizing 0(0;^) 
corrections are colour suppressed and are thus expected to be small as e. g. explicitly 
shown for B — > D^*^Tr~ in As concerns the nonperturbative effects we also write 
down the L/T composition of the nonperturbative 0(l/m^) contribution to the D* 
rate as well as to the spin Ds rate using results of |[T0[| . 

At the Born term level Aleksan et al. found ri/rj' = 1.823 (using their mass values 
mft = 4.85 GeV and mc = 1.45 GeV). It would be interesting to see how radiative correc- 
tions and the nonperturbative contributions affect this ratio. In the corresponding case 
t—>b + (with mf, = 0) we found earlier that the ratio Vl/Vt is shifted downward 
by the 0{as) radiative corrections by an amount of 3.5% ||Tl[| . 

2 Angular decay distributions 

The longitudinal and transverse content of the diagonal density matrix of the Dl~ 
(or its charge conjugate state Dl'^) can be determined by analysing the angular decay 
distribution of its subsequent decay into D*~ — > D~ + 7 and Dl~ D~ + 7r°. The 
branching ratios into these two principal channels are given by (94.2 ± 2.5)% and 

"'^The 0{as) corrections to the spin 1 piece of the weak current keeping both quark masses finite 
had been calculated before in [^|-^. The 0{as) corrections to the spin piece of the weak current 
can be deduced from the corresponding calculation for t ^ b + Q . Latter result had also been 
used in a calculation of the 0{as) radiative corrections to h ^ c + + Vr where the spin piece 
enters because the r-mass cannot be neglected in this process. 
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(5.8 ± 2.5)% [0, respectively. In terms of the diagonal density matrix elements pn 
{■m = 0{L), ±1(T)) of the D*~ the polar angle distribution is given by 



The hm are the decay amplitudes of the decays D* + •y (m = ±l) and D* —>■ 

D~ + -7r° (m = 0) where the m are the magnetic quantum numbers of the D*~ in the 
decay frame. The d^^^i{0) are the usual Wigner c?-function and 6 is the polar angle 
of the D~ in the D*^ rest frame (measured with regard to the original momentum 
direction of the D*~) as shown in Fig. 1. One thus obtains the polar angle decay 
distributions 

d^B^X,+Dr{^DT+-/) _T3T3/n*- . n- , ^^^^^^ , ^„r.2n\-r , ^ .^^2/ 



dcosO BR(Z)r^Z);+,)^-(l+oos^O)r, + jsm^Or,j (2) 

and 

^-^ = BR{D, ^D, +vr)^-sm^r, + -cos^r,J. (3) 

Considering the fact that the upcoming 5-factories will be producing upward of 
lOK BB pairs per day and that the inclusive branching ratio of the -B's into -D*'s is 
expected to lie around 0(5%) it should not be too difficult to experimentally determine 
the angular coefficients of the two decay distributions and thereby the L/T content of 
the D*. 

3 Born term rates and 0{cXs) radiative corrections 

Let us begin by writing down the Born term level results for b—>- {D^,D*^) + c (see 
Fig. 2a). We shall closely follow the notation of Aleksan et al. ^ throughout. For easy 
comparison with the numerical results of we shall also adhere to their numerical 
parameter values. One has 



Svr mi \ [mi — m^^Y / 

r,^^(.^D. +0 - ^IHcKJ fn-t ^2 (,1+ (^2 _ ^2)2 j^^^J^i' (5) 

r. i^-r>t +0 = — iHcKJ fn, (1- (^2_^2)2 (6) 



In Eqs. (§-|D fos and foi denote the pseudoscalar and vector meson coupling con- 
stants defined by (L)7|A'^|b) =i/z),p^^ and {D*-\V^'\Q) = foimo^e*^ , respectively. The 
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Kobayashi-Maskawa matrix element is denoted by V^^qj ^^e po^ and po* are the 
three-momenta of the Ds and D* in the b rest system. The parameter ai is related to 
the Wilson coefficients of the renormalized current-current interaction and is obtained 
from a combined ffi of several decay modes (|ai| = 1.00 ± 0.06) [Q. Note that the 
structural similarity of the rate formulae for the decay into Dg and the longitudinal D* 
is an accident of the Born term calculation and does not persist e.g. at higher orders 
of as- 

Using /d. =230 MeV and /d* =280 MeV as in H, tb = 1.6 ps, He = 0.04, K. = 0.974 
and the central value for oi one arrives at 

BR,_^-+, = 3.2% BR,_,^.-+, ^ 6.8%. (7) 

Summing up the Ds and D* modes one arrives at a branching ratio of 10% which is 
consistent with the measured value BR{B ^ DfX) = (10.0 ±2.5)% |T^ if one assumes 
that the above two rates saturate the inclusive rate into Df. 

Next we turn to the 0{as) radiative corrections. As explained in 0], the radiative 
gluon corrections connect only to the b and c legs of the parton decay process b — > 
{Ds, -D*) -|-c because of the conservation of colour (see Fig. 2b, 2c and 2d). As remarked 
on earlier the radiative corrections for the spin 1 piece are then identical to the radiative 
corrections calculated in [^j-^ or in [l^ where the process t —>-W^ + b was considered 
keeping rUh^O. In [|T^ we separately computed the radiative corrections to longitudinal 
(L) and transverse (T) W^+'s in the decay process and thus these results can directly 
be transcribed to the present case. 

The two L and T pieces can be projected from the hadron tensor by use of the 
projection operators^ fl^ (we use T = (L + T) — L since (L + T) is simple compared 
to either L or T) 



= (- + 



The scalar spin piece can be obtained with the projector 

Fr = (^). (10) 

The four-momentum of either the Dg or the D^ is denoted by g^. The magnitude 
of the three-momentum of is given by |g| = \/ (f — q^- We shall not dwell much on 



the details of our calculation in this short communication but refer to ||13| for technical 
details. Let it be said that we use a gluon mass regulator to regularize the infrared 
singularities differing from Aleksan et al. who use dimensional regularisation instead. 



^Here we use the notation "T" ("transverse") rather than the notation "[/" ( "unpolarized trans- 
verse") used in jl^. 
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We shall present our 0(0;^) results in a form where the respective Born terms vf^ are 
factored out from the 0{as) result. Including the Born term and the nonperturbative 
0{l/ml) contributions to be discussed in Sec. 4 we write with f 5 :=r5'/r|''*, tL^T,L-= 
i L+r,L/i and i _5 



r, = r,^'^^(l + Cf—T, + K,ai + G,h), (11) 



TT 



where i = S, L + T, L. Kb and Gh are the expectation values of the kinetic energy and 
the chromomagnetic interaction of the heavy quark in the B meson, respectively. 

To begin we list the reduced 0{as) rates Fj. For the reduced scalar spin rate F^ 
we obtain 

= 2E + 1 - |^lnK^,) + A-V25/{|^i^4^1^(^)+^(l-^)' X 
(l-x+y)(l+a; + y)lnKw^)-^^/2(l-x2)(l+y2)7^(_2_l) + l(l+^2) ^ (^3) 

where we have defined a Born term-like scalar rate by 

Bs = {l-y'f-x\lW). (13) 

The other variables and functions appearing in Eq. (|12[) are explained at the end of 
this section. For the total reduced spin 1 rate we obtain 

1 ( 1 (^ _,,2^ \3/2 1 

W = 2S - --^\n{w,w,) + \-'"B,l^[-^^^ ln(y) + _(A + 6x2y) x 

{l-x+y){l+x+y) ln(wlw;^)-i^/2(l_a;2^(^^2x2+^/2)7^(_2,-l) + ^((l+a;^- 
2x^) + (4-3x2)l/2 + 3y4)7^(_,__^-l(3-2x2 + 3y2)7^(o,-l)-^ (1 + 2x^+1/2) x 

l/'5(o,o) + ^(l + 2x2+|/2)5(i,o)}, (14) 

with 

BL+T = {l-y^?+x\l-2x^+y^). (15) 
Finally, the longitudinal piece of the reduced spin 1 rate is given by 

f L = 2 S - 1^ Hw^w,) + \-'''B,^{\ ^^-y'l^"\ n{y) + y {l-yf x 
4 x^ 1 2 x"^ 

(l-X + y)(l+X + y) HwiW^) - i(l-x2)3(l+y2) y2 7^(_2^^) + i(l-X=') X 

5 



((l-a;2)^ + (6+x2-3x^)?/2+(5 + 3x2)?/^)7^(^l,l)-^((5-2x2-7x^+4x^) + 
(12-33x2 + x^)l/2 + (7+x2)l/^)7^(o,l) + ^((7-31x2+4x^) + (10+x2) + 
2,y^)n^i,i) - ^(l+y2)7^(2,l) - i((H-10x=^-llx^) + (l+x2)V)y''5(o,2) + 
^((l + 10x2-llx^) + (3-4x2+x^) 1/2 + 2 (l+x2)y^)5(i,2)-^((2-6x2) + 

„2\ „.2 I „.4 



(3+2x2) l/2+y^)5(2,2)+2(l+i/2)5(3 



2) f, 



where 



BL = Bs = {l-yy-x\l+y 



(16) 
(17) 



The contribution denoted by S is the finite remainder of the Born term type one- 
loop contribution plus the soft gluon contribution. It is given by 



l-x^+y'^ 
AV2 



Wi 



Li2 1 1 ) - Li2 ( 1 - ) - Li2 ( 1 - wi w^) + - In 



w'l 



In 



xy^ 



— In Wi In 



i-i.,/ 

2 





x^~2yi 




12 ' Ai/2 J 



1-wf 

2 



4A1/2 '^U„ 



(1^ 



where we use the abbreviations 

{l-x^+y^-X^/^)x 
(l+x2-|/2 + AV2)y' 



Wi 



{l-x^+y^-\'/^)x 

(l + x2-|/2_Al/2)2/' 



(19) 



and x = m^{t)/mi, and y = mc/mb- The kinematical factor A is defined by A = l + x^ + 



y'^ — 2 x'^ — 2 y'^ — 2 x'^y'^ such that p^(*) = A^/2_ 

The reduced 0(0;^) rates are given in terms of a set of tree graph phase space 
integrals 7l(m,n) and which are defined by 



7e 



(m,n) 



n/2 



A 



(m,n) 



m 



A 



n/2 



1— x^+z— A 



1/2 



(20) 



where Xz = l+x'^+z'^—2x'^—2z—2zx'^. Their solution can be obtained using techniques 
similiar to the ones discussed in 113. With the three abbreviations 



Afi := Li2(ux)-Li2( 



M2 ■= -lnMln(l+x)+ln 



u—x 



In 



(m — l)(l+x)/ \m(1 — Mx) 



u—x 



+ 



(21) 
(22) 



6 



- Ll2 - +L12 r +L12 



u/ \u{u—x)J \ u — x 

Ms -ln.ln(l-.)-ln ((!l±iKlZ:^) In f-^- ) + 

V u — x / \u [l—ux) ' 

- L12 -- +L12 ^ ^ +L12' 



u/ V u—x / \ u—x 
one has 

^5(0,2) = -^(M-M), 

•^(i-s) 2x •^3+.A/i, 

5(2,2) = -^^^AA2 + ^^^AA3+2(l+x2)M + A^/^-2x2ln« + 
^ ^ 2x 2x 

\U[1 — UX)/ 

5(3,2) = -^^^A/'2+^^^A/'3+(3+x^)(l+3x^)M + ^(9 + 13x^ + y^)AV2+ 
- (6+5 x^) x' lnu + ^ {4{l+x') -y') In ( J . 

The non-logarithmic integrals are given by 
7^(_2-l) = — A^/^+lnw-^ ^In 



1 — Vl— -ua;/' 
7e(_i,_i) = -A^/2-(l+a:2)ln«+(l-x2)ln(-^), 

7^( o,-i) = \{l+x''-y^)X'/^-2xHnu, 

_ 1 rA^/^ f u-x 

^M. ^) = rr^^^(ir^)' ^(0,1) = 1^^' n^,,^^=-X'/' + {l+x')lnu, 
7^( 2, 1) = -^(3+3x'+y')A^/' + (l+4x'+x^)lnM, 

l+x^-y^+VA 

where . 

2x 



4 Nonpertubative contributions 



When one uses the operator product expansion in HQET one can determine the nonper- 
tubative corrections to the leading partonic 6— >c rate. The nonpertubative corrections 
set in at 0(l/m^) and arise from the kinetic energy and the chromomagnetic interaction 
of the heavy quark in the heavy hadron. The strength of the kinetic and chromomag- 
netic interactions are parametrized by the expectation values of the relevant operators 
in the B system and are denoted by Kb and Gb, respectively. The nonpertubative 
contributions to the spin and spin 1 rates including the L / T separation have been 
calculated in and can be taken from there. One has 



where 



S : as = -1, 

bs = {BsX)-'[-{l-yy{l-5y^)+x\3-7y'-lly^+15y' 

-x^{7+l0y^+ 15?/^) +5x^(1+7/2) 
16 

L: aL = -l--x^Bl\ 

bL = (35iA)-^[-3(l-y')3(l-5y')-a;2(7-27y' + 65y^-45i/^ 
+x^(27+3%2-45?/^)-a;6(17-15y2) 

T : ar = -1 + —x^B^^, 
3 

br = 2x^(35tA)"^ {l-y^){5-4:y^ + l5y^)+x^{9 + l0y^+A5y^ 
-x^{29+A5y^) + l5x^ 



Bt = Bl+t -Bl = 2x\l-x^ + y^). 
For our numerical evaluations we use Kb = 0.013 and Gb = —0.0065 as in [TO]. 



(36) 
(37) 

(38) 
(39) 

(40) 
(41) 



(42) 



5 Numerical results 

Using mft = 4.85 GeV, mc = 1.45 GeV, rriD^ 
«s(^fe) = 0.2 we obtain for b^c 



1968.5 MeV and m^* =2112.4 MeV and 



Ts = (1 - 0.09638- 0.013 +0.00467), 

Tl = 0.6459(1 - 0.1103 - 0.03413 + 0.00966), 
fr = 0.3541 (1 - 0.1079 +0.02553- 0.02769), 
fi+r= (1 - 0.1095 - 0.00859- 0.01803). 



(43) 
(44) 
(45) 
(46) 
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The radiative corrections reduce the rates by about 10%, where the reduction is rather 
uniform for the four different rates. The nonperturbative corrections range from 0.5% 
for the chromomagnetic correction to Ts to a maximal 3.4% for the kinetic energy 
correction to Tl with no uniform pattern in their contributions. At the Born term 
level the transverse/longitudinal composition is given by rj-Zr/, = 0.55. This ratio is 
shifted upward by the insignificant amount of 0.3% through the radiative corrections. 
Adding all corrections one finds a 3% reduction in the ratio. 
For the b^u transitions with mu = 0, i.e. y = we have 

Ts = (1 -0.1694- 0.013 +0.00751), (47) 

Tl = 0.7250(1 -0.1777- 0.02923 + 0.01414), (48) 

fr = 0.2750(1 -0.1150 + 0.02978- 0.02348), (49) 

Tl+t= (1 -0.1605 + 0.00055- 0.00934). (50) 



In the h^u case the dominance of the longitudinal rate is more pronounced. At the 
Born term level one finds r-r/ri = 2x^ = 0.38. The radiative corrections are no longer 
as uniform as in the h^c case. Whereas the radiative corrections to Ts-, Tl and r^+T 
amount to 16% — 17%, the radiative correction to the transverse rate Tt is only 11.5%. 
Thus the ratio Tt/Tl is shifted upward by 7.6% by the radiative corrections. Adding 
up all corrections one finds a 10.4% upward shift for this ratio. Let us mention that 
our 0{as) results on F^+t and Ts numerically agree with the results of [0] for both the 
h^c and h-^u transitions. 

As emphasized in the introduction the conclusions drawn in this paper on the ra- 
diative corrections are tentative in as much as there are also nonfactorizing 0{as) 
contributions which have not been included in our analysis. Although the nonfactoriz- 
ing 0{as) contributions are colour suppressed and thus expected to be small it would 
nevertheless be worthwhile to try and estimate the nonfactorizing 0{as) contributions 
along the lines of and 0. 

The last point we want to discuss are the inclusive decays B—^Xc+ {ir^ , p^) which 
can also be induced by the diagrams Fig. 2 when the c^s transition in the upper leg 
is replaced by a -u^c/ transition. Using = 132 MeV, fp- =216 MeV and V^rf = 0.975 
one finds the Born term branching fractions BRh^„-+c — 1-6% and BRb_^p-^c — 4.6%. 
In the latter case the rate is dominated by the longitudinal contribution since q^ = 7ii?p 
is not far from = where the rate would be entirely longitudinal. In fact one finds 
Tt/I^l = 0.067. It is important to note that the diagrams Fig. 2 are not the only 
mechanisms that contribute to the inclusive decays B^Xc+ {tt~,p~). Additional vr" 
and p~ mesons can also be produced by fragmentation of the c-quark at the lower leg.Q 
As concerns the p~ mesons resulting from the fragmentation process they would not 

concerns the inclusive decays B Xc + (D^ , D* ) the possibiUty of producing extra and 
D*~ mesons through fragmentation of the c-quark is ruled out for kinematic reasons. 
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be polarized along their direction of flight. This lack of polarization as compared to 
the strong polarization of the p mesons from the weak vertex could possibly be used 
to separate p~ mesons coming from the two respective sources. 
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Figure 1 : Definition of polar angle 9 in the inclusive decay 5 — > Xc+-D* (— > D~ +7 or 7r°) . 
The polar angle 9 is defined in the D*" rest frame relative to the direction of the D*~ 
in the B rest frame. 



Figure 1 




Figure 2: Leading order Born term contribution (a) and 0{as) contributions (b,c,d) to 
b^c+{D7,D*-). 



Figure 2 
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